OPEN 3 ACCESS Freely available online 



•0-PLOS I ONE 



A Model of Two-Way Selection System for Human 
Behavior 

Bin Zhou 1,2 , Shujia Qin 3 , Xiao-Pu Han 4 , Zhe He 1 , Jia-Rong Xie 1 , Bing-Hong Wang 1 ' 5 * 

1 Department of Modern Physics and Nonlinear Science Center, University of Science and Technology of China, Hefei, China, 2 Suzhou Institute of Technology, Jiangsu 
University of Science and Technology, Suzhou, China, 3 State Key Lab of Robotics, Shenyang Institute of Automation, Chinese Academy of Sciences, Shenyang, China, 
4 Institute of Information Economy and Alibaba Business College, Hangzhou Normal University, Hangzhou, China, 5 College of Physics and Electronic Information 
Engineering, Wenzhou University, Wenzhou, China 



Abstract 

Two-way selection is a common phenomenon in nature and society. It appears in the processes like choosing a mate 
between men and women, making contracts between job hunters and recruiters, and trading between buyers and sellers. In 
this paper, we propose a model of two-way selection system, and present its analytical solution for the expectation of 
successful matching total and the regular pattern that the matching rate trends toward an inverse proportion to either the 
ratio between the two sides or the ratio of the state total to the smaller group's people number. The proposed model is 
verified by empirical data of the matchmaking fairs. Results indicate that the model well predicts this typical real-world two- 
way selection behavior to the bounded error extent, thus it is helpful for understanding the dynamics mechanism of the 
real-world two-way selection system. 
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Introduction 

Human-initiated systems always run in a complex way. In the 
past ten years, related work mainly focused on the temporal and 
spatial distribution characteristics of human activity patterns. 
Because of the complexity of human behavior, many underlying 
mechanisms have not been discovered yet. The two-way selection 
scenario among humans is one of the complicated but common 
phenomena in daily life. It happens in the processes like choosing a 
mate between men and women, making contracts between job 
hunters and recruiters, and trading between buyers and sellers. In 
a sense, two-way selections can be regarded as the base of building 
many social relationships. Generally, the participants in a two-way 
selection process are first classified into two groups by their natural 
status. Then they observe, study the factors of the people on the 
other side, and finally make their choices. For instance, in the case 
of marriages, one's appearance, personality, wealth, and sense of 
humor, are prevalently taken into consideration. Besides the 
individual characters, impersonal factors also exert an influence, 
e.g. the member totals on each side and their ratio. How many 
characters will be inspected and chosen deeply affects the result of 
a selection process. However, usually it is difficult to compare and 
to distinguish these characteristics quantitatively even qualitatively 
through traditional methods, such as psychological tests and social 
surveys. 

The well-known marriage game in statistical physics has been 
researched in these papers [1-6], whose main novel concept is the 
stability of marriages. This view point aims to find a stable 



matching between the two sets of men and women. Such a model 
results in the destiny that every one in the sets gets married and the 
final marriage relationships are "stable". However, the internal 
mechanism of a two-way selection system can be modeled in 
another way: not all of the participants have to get married in one 
trial of the processes, i.e. some of them would be successful in 
matching but the others not. This mechanism would render 
assistance to some social problems, such as the prediction of the 
total of friendships or other gregarious relations [7-15]. In this 
paper, we present a model for two-way selections to investigate the 
factors influencing the matching rate. The data of matchmaking 
fairs are analyzed to support our model. Based on this model, the 
method of estimating the number of factors impacting people's 
decisions is also proposed. 

The Model and Analytical Results 

Our model of the two-way selection is stated as follows: 

i) The system has two sets of agents, A and B, respectively 
amounting to k\ and lc2. 

ii) The zth agent in set A (or set B) has its own character 
denoted by Ca, (or eg,). Correspondingly, the character the 
zth agent attempts to select is denoted by Sa, (or Sg,). 

iii) The agents' characters are denoted by integers without loss 
of generality. Assume the characters has n types, i.e. Ca,, Cb j? 
s a „ s Bj eQ, n = { 1,2, ...,«}, fe{l,2,...,ifci}, 
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je{l,2, . . . ,k 2 }. In one trial of the model, ca,, cb p sa„ sbj 
pick an element in S following the uniform distribution, 
iv) The condition of successful matching of two agents Aj and 
Bj is ca,=sbj and sa, = cb s - That is, when agent A,'s 
character meets agent B/s requirement and vice versa, agent 
At and agent Bj have a successful matching. 

For given k\, k 2 and n, the expectation E of the total number of 
matching pairs in the model is 



E(k\,k2,n) 



i = 0 \J=i+\ 



E C 'ki-l CJ k 2 [- 



1 \ i+j (n 2 - l\ kl +k 2~'~ J ~~ 



7 = 0 



n 2 -l 

2 



k\ + k 2 —i—j — 1 



According to 



lim C k x l (\-x) k -''= e -^, 

00 x' 

^^--=BesselI(0,2^), 

00 x ! 

E ni =\/^BesselI(l,2^), 



Because in this case k 2 is very large, Equation (5) can be further 
simplified as 



Define 



Exk x . 



= ki t= r? p _ 2E 
^ k\ ' k\ ' fci +&2 ' 



(6) 



(7) 



where r] denotes the ratio of ki Xok\\ t, denotes the ratio of « 2 to 
k\ ; P denotes the estimated ratio of successful matching pairs to 
the average number of two type agents. Then Equation (2) can be 
transformed into 



1 _ e ~H (Bessell (fJ, ^ + J- Bessell ^1 , ^ 



Equation (3) can be written as: 
2r\ 



(8) 



(9) 



Equation (6) can be written as: 



where X = xk, Bessell^^^/x) and Bessell(l,2 % /x) are the 
modified Bessel functions of the first kind, the expectation in (1) 
can approximate to 



Exki [l-e~ x i - x 2 (Bessell (o,2\/khj 



+ W --Bessell ( 1,2 y Xih 



(2) 



where X\ =k\/n 2 , ).2=k2/n 2 . 

Due to the symmetry of k\ and k 2 in (1), without loss of 
generality, we just study the k\ < k 2 case under three conditions: 
ki<k 2 «n 2 , k\<n 2 «k 2 , and n 2 «k\<k 2 . When k\<k 2 «n 2 , 
resulting in k\ ->0,A2->0, calculating the zeroth power term and 
the first power term of (2) obtains 



Ek \k\k 2 . 



When k\ <n 2 «k 2 or n 2 «k\ <k 2 , according to 



lim Bessell(0,2x) = lim Bessell(l,2x) = 



Equation (2) can be simplified as 



2y / 7Tx' 



(3) 



(4) 



Exki—ki 



,-(\/h-^r 2 ) 2 

2^JtQ-XiX 2 




(5) 



f7+l 



, 1 <c;«ii or ^«\<r\. 



(10) 



Figure 1 shows the comparison between the analytical 
predictions of (1) and the simulation results. Figure 2(a) shows 
the comparison between the analytical predictions of (9) and the 
simulation results under the condition 1 < r] « £, and displays a 
power-law relation with the exponent —1 between P and £. 
Figure 2(b) shows the comparison between the analytical 
predictions of (10) and the simulation results under the condition 
1 < £ « rj and displays a power-law relation with the exponent — 1 
between P and rj; Figure 2(d) shows the comparison between the 
analytical predictions of (10) and the simulation results under the 
condition £ « 1 < r] and displays the same power-law relation to 
the result in Figure 2(b). The above analytical predictions and 
simulation results are consistent with each other. That is to say all 
analytical results are reliable. 

Consider a special case k\ =k 2 =k, resulting in rj=l. On the 
one hand, Equation (9) can be simplified as 



(11) 



The relation between P and £ approximates a power law with the 
exponent — 1, and this case is shown in Figure 2(c). Equation (10) 
can be simplified as Px 1, suggesting that almost all of the agents 
can match successfully under the condition £ « 1 < r\. On the other 
hand, because the condition k\=k 2 = k results in X\ =X 2 , from (5) 
we can obtain 
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Figure 1. The comparison between the analytical predictions and the simulation results. In the two sub-figures, the parameter k\ = 100; r\ 
is assigned to values from 1 to 1000; £ is assigned to values from 0.1 to 1000. (a) shows analytical predictions of (1); (b) shows the simulation results. 
doi:10.1371/journal.pone.0081424.g001 



E&k-J—Vk. (12) 

V 71 

The second term of ( 1 2) is the number of the agents that can not 
successfully match in type A or type B. The larger k is, the smaller 
\J n 2 1 n\fk /k is. In reality, this is a result of fluctuation. The total 
combinations of the "own" state and the "expecting" state for an 
agent have n 2 possibilities in the model. In theory, the expected 
times of each state appearance is k/n 2 . However, due to the 
fluctuations, almost all frequencies of every state appearance 
deviate around k/n 2 . As a result, some agents can not successfully 
match. The number of times that each state may appear obeys the 
binomial distribution. The fluctuation is closely related to the 



standard deviation, according to the binomial theorem and 
standard deviation formula, we can obtain that the standard 
deviation equals yjk(\ — 1 /n 2 )/n 2 , which is directly proportional 
to \fk. Thus, the number of the agents that can not successfully 
match is also proportional to \fk. It explains the relationship 
between the second item of (12) and \fk. From (12), we know that 
the proportionality coefficient is \J n 2 Jn. 

The Verification Between the Model and 
Experimental Data 

As the mate choosing between men and women is a typical real- 
world two-way selection system, eighty-two reported records of 
matchmaking fairs are analyzed to verify our model. Due to the 
uncertainty of approximation in these reports, we classify the data 




Figure 2. The comparison between analytical predictions and the simulation results in the log-log plots. In the four sub-figures, the 
parameter k\ = 100; the squares are the simulation data; the solid lines are analytical predictions. In (a), (j = 2; { is assigned to values from 100 to 1000; 
the solid line is obtained from (9). In (b), {=2; is assigned to values from 10 to 1000; the solid line is obtained from (10). In (c), r\= 1; { is assigned to 
values from 10 to 100; the solid line is obtained from (1 1). In (d), c = 0.1; 17 is assigned to values from 10 to 1000; the solid line is obtained from (10). 
doi:10.1371/journal.pone.0081424.g002 
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Table 1. The data of matchmaking fairs. 



Website no. 


Original descriptions 




Total participants K 


Matched pairs E 


Matching ratio P 




joined 


matched 








01 


13 


3 


13 


3 


0.4615 


02 


18 


6 


18 


6 


0.6667 


03 


20 


1 


20 


1 


0.1000 


04 


21 


6 


21 


6 


0.5714 


05 


25 


5 


25 


5 


0.4000 


06 


26 


3 


26 


3 


0.2308 


07 


26 


4 


26 


4 


0.3077 


08 


30 


2 


30 


2 


0.1333 


09 


32 


4 


32 


4 


0.2500 


10 


36 


6 


36 


6 


0.3333 


1 1 


36 


6 


36 


6 


0.3333 


12 


38 


6 


38 


6 


0.3158 


13 


40 


8 


40 


8 


0.4000 


14 


>40 


3 


42 + 2 


3 


0.1432 ±0.0068 


15 


>50 


3 


53 + 3 


3 


0.1136±0.0064 


16 


=60 


5 


60 + 3 


5 


0.1671 ±0.0084 


17 


>60 


4 


63 + 3 


4 


0.1273 ±0.0061 


18 


>60 


5 


63 + 3 


5 


0.1591 ±0.0076 


19 


>60 


— 10 


63 + 3 


10±1 


0.3197±0.0470 


20 


72 


11 


72 


1 1 


0.3056 


21 


80 


5 


80 


5 


0.1250 


22 


80 


10 


80 


10 


0.2500 


23 


80 


13 


80 


13 


0.3250 


24 


80 


18 


80 


18 


0.4500 


25 


-100 


5 


95±5 


5 


0.1056±0.0056 


26 


99 


8 


99 


8 


0.1616 


27 


100 


5 


100 


5 


0.1000 


28 


~100 


7 


100±5 


7 


0.1404 ±0.0070 


29 


>100 


16 


110±10 


16 


0.2933 ±0.0267 


30 


150 


7 


150 


7 


0.0933 


31 


-200 


8 


190± 10 


8 


0.0844 ±0.0044 


32 


-200 


22 


190± 10 


22 


0.2322±0.0122 


33 


-200 


4 


200± 10 


4 


0.0401 ±0.0020 


34 


206 


10 


206 


10 


0.0971 


35 


>200 


7 


210 + 10 


7 


0.0668 ±0.0032 


36 


>200 


8 


210+ 10 


8 


0.0764 ±0.0036 


37 


>200 


38 


210+ 10 


38 


0.3627±0.0173 


38 


216 


19 


216 


19 


0.1759 


39 


>240 


22 


252 + 12 


22 


O.1750±0.0083 


40 


-258 


>10 


258 + 13 


11±1 


0.0859±0.0121 


41 


-300 


4 


285 + 15 


4 


0.0282±0.0014 


42 


=300 


8 


300 ± 1 5 


8 


0.0535±0.0027 


43 


>300 


>10 


315 + 15 


11±1 


0.0703 ±0.0097 


44 


>300 


32 


315+15 


32 


0.2036 ±0.0097 


45 


400 


-20 


400 


19±1 


0.0950 ±0.0050 


46 


>500 


3 


525 ±25 


3 


0.0115±0.0005 


47 


>500 


8 


525 ±25 


8 


0.0306±0.0015 


48 


>500 


>10 


525 ±25 


11±1 


0.0422 ±0.0058 
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Table 1. Cont. 





Website no. 


Original descriptions 




Total participants K 


Matched pairs E 


Matching ratio P 




joined 


matched 








49 


-600 


-40 


570 + 30 


38±2 


0.1341 ±0.0141 


50 


>600 


>78 


630 ±30 


78 


0.2482±0.0118 


51 


-800 


58 


760 ±40 


58 


0.1531 ±0.0081 


52 


>800 


>20 


840 ±40 


21 ±1 


0.0502 ±0.0047 


53 


-1000 


=20 


950 ±50 


20±1 


0.0423 ±0.0043 


54 


-1000 


58 


950 ±50 


58 


0.1224±0.0064 


55 


-1000 


64 


950 ±50 


64 


0.1351 ±0.0071 


56 


"1000 


12 


1000 ±50 


12 


0.0241 ±0.0012 


57 


-1000 


15 


1000 ±50 


15 


0.0301 ±0.0015 


58 


"1000 


-100 


1000 ±50 


95±5 


0.1910±0.0196 


59 


>1000 


3 


1050 ±50 


3 


0.0057 ±0.0003 


60 


>1000 


4 


1050 ±50 


4 


0.0076 ±0.0004 


61 


>1500 


48 


1575 ±75 


48 


0.0611 ±0.0029 


62 


>1500 


>100 


1575 ±75 


105±5 


0.1339 + 0.0127 


63 


>1600 


31 


1680 ±80 


31 


0.0370±0.0018 


64 


>2000 


-100 


2100±100 


95±5 


0.0909±0.0091 


65 


>2000 


>113 


2100±100 


119±6 


0.1139±0.0111 


66 


-3000 


-100 


2850 ±150 


95±5 


0.0670 ±0.0070 


67 


=3000 


186 


3000 ±150 


186 


0.1243 ±0.0062 


68 


>3000 


>200 


3150± 150 


210±10 


0.1339±0.0127 


69 


>4000 


>500 


4200 ±200 


525±25 


0.2511 ±0.0239 


70 


-5000 


108 


4750 ±250 


108 


0.0456 ±0.0024 


71 


>5000 


218 


5250 ±250 


218 


0.0832 ±0.0040 


72 


>5000 


231 


5250 ±250 


231 


0.0882 ±0.0042 




>5000 


237 


5250 ±250 


237 


0.0905 ±0.0043 


74 


>6000 


>270 


6300 ±300 


284±14 


0.0906 ±0.0088 


75 


— 1 0000 


28 


9500 ± 5U0 


28 


0.0059 ±0.0003 


76 


-10000 


= 100 


9500 ±500 


100±5 


0.0212±0.0022 


77 


-10000 


-2000 


9500 ±500 


1900±100 


0.4022 ±0.0422 


78 


> 10000 


=400 


10500 ±500 


400 ±20 


0.0766 ±0.0075 


79 


> 10000 


>1000 


10500 ±500 


1050±50 


0.2009±0.0191 


80 


=16000 


>700 


16000 ±800 


735 ±35 


0.0923 ±0.0090 


81 


> 16000 


>600 


16800 ±800 


630 ±30 


0.0753 ±0.0072 


82 


>50000 


>3000 


52500 ±2500 


3150± 150 


0.1206±0.0115 















Note: = denotes "about", — denotes "nearly", and > denotes "above". 
doi:1 0.1 371 /journal.pone.0081 424.t001 



into three categories with specified possible ranges according to 
their descriptions: i) "nearly x" (possible range 0.95x + 0.05x); ii) 
"about x" (possible range 1.00x + 0.05x); iii) "over x" (possible 
range 1.05x + 0.05.\:). The full list of the data records is shown in 
Table 1. All data of matchmaking fairs are collected from the 
websites shown in Table 2. 

In our model, n is an internal parameter needed to be measured. 
Because a news report (descried as an experiment below) generally 
includes only the total of participants and the number of successful 
matching pairs, the male-female or female-male ratio r\ defined in 
(7) should be estimated first. Under the assumption k\<k2, the 
lower bound of rj is >7 min = 1 , and once the total of participants 
K = k\ + ki and the number of matching couples E is determined, 
the upper bound of rj in that experiment is known: r]<(K — E)/E. 



Let jV be the number of experiments, f/ max> , be the upper bound of 
v\ in the ith experiment, and H = {i? max-1 ,>7 max> 2, ■ ■ • ,'7max,jv} be the 
set of all upper bounds. By processing N = 82 experiments in 
Table 1, we obtain min(7/) = 2 and max(//) = 365.67. Consider 
the least square criterion for fitting the model and the experimen- 
tal data 

N 

Q(n, tl )=J2(P e ,i-PrAn,r 1 )) 2 (13) 
1 = 1 

where P e ,i denotes the experimental data in the zth experiment 
and P,j(n,ii) denotes the corresponding theoretical value calcu- 
lated by (8), and the reality that in a matchmaking fair the 
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Table 2. Data 


sources of matchmaking fairs. 




Website no. 


Website name of matchmaking fairs 


01 


http://www.cdb.org.cn/newsview.php7id = 6359 


02 


http://sd.people.com.en/n/201 2/0827/cl 8371 8-1 7407663.html 


03 


http://news.carnoc.eom/list/1 83/1 83765.html 


04 


http://bbs.tiexue.net/post2_5756757_1.html 


05 


http://www.shxb.net/html/20110516/20110516_278862.shtml 


06 


http://news.qq.eom/a/201 0051 1/000472.htm 


07 


http://bbs.ganxianw.com/thread-46316-1-1.html 


08 


http://www.wzrb.com.cn/article321273show.html 


09 


http://www.wccdaily.com.cn/epaper/hxdsb/html/2012-05/14/content_448239.htm 


10 


http://wbnews.sxrb.com/news/ty/1372966.html 


11 


http://www.nbmz.gov.cn/view.aspx7id = 16595&AspxAutoDetectCookieSupport= 1 


12 


http://nb.people.com.en/GB/200892/1 6491 824.html 


13 


http://cq.cqwb.com.cn/NewsFiles/201203/25/921397.shtml 


14 


http://www.sc.chinanews.com.cn/my/data/html/201212/32619.html 


15 


http://www.1 6466.com/info_detail.htm7id = 36526 


16 


http://www.ncnews.com.cn/ncxw/shxw/t20121 1 12_9431 14.htm 


17 


http://news.wzsee.com/2012/0502/130061.html 


18 


http://news.hexun.com/2012-08-27/145162214.html 


19 


http://www.douban.com/group/topic/28145139/ 


20 


http://zhuanti.1 0yan.com/zt/other/sdcms/html/xqj201 2/xiangqindongtai/1 377.html 


21 


http://www.wlmqwb.com/3229/syzt/hdzt/seven/201007/t20100719_1287834.shtml 


22 


http://www.dpcm.cn/html/news/shehui/20121211/8a485b96f289e038.htm 


23 


http://www.zhaogejia.com/News/Show/166 


24 


http://cq.cqnews.net/shxw/shwx/200909/t20090928_3635782.htm 


25 


http://a.jiaodong.net/jiaoyou/detail/?/201 2071 71 3471 5.htm 


26 


http://www.dllake.com/testurl/news/news.asp7id = 1874 


27 


http://fj.qq.eom/a/201 2041 3/000073.htm?pgv_ref = aio201 2&ptlang = 2052 


28 


http://epaper.lnd.com.cn/html/bdcb/201 1011 8/bdcb635730.html 


29 


http://news.zh853.com/NewsShow-22166.html 


30 


http://news.163.com/11/1123/08/7JHJ4PP100014AED.html 


31 


http://news.xinmin.cn/rollnews/201 1/05/03/1 0539635.html 


32 


http://www.0523qq.com/forum.php7mod = viewthread&tid = 2779 


33 


http://news.ycw.gov.cn/html/2012-04/28/content_15150376.htm 


34 


http://epaper.lnd.com.cn/html/bdcb/201 1011 8/bdcb635730.html 


35 


http://www.cqwb.com.cn/NewsFiles/201 005/30/201 0293006291 035471 6.shtml 


36 


http://bddsb.bandao.cn/data/20120827/html/53/content_2.html 


37 


http://www.zhaogejia.com/News/Show/150 


38 


http://3g.3xgd.com/news/play.asp7NewslD = 80975 


39 


http://wed.cnhan.com/hjb/201 2-1 2-03/3900.html 


40 


http://xt.fangyuan365.com/article/List.asp7ID = 8708 


41 


http://cq.cqnews.net/shxw/201 2-1 1/1 2/content_214321 70.htm 


42 


http://xt.fangyuan365.com/article/List.asp7ID = 11694 


43 


http://roll.sohu.com/20120625/n346389451.shtml 


44 


http://www.ijxjj.com/article/article_12773.html 


45 


http://news.xinmin.cn/shehui/201 3/02/1 6/1 8638248_2.html 


46 


http://www.sz120.com/xwdt/ynxw/22205/ 


47 


http://www.sc.xinhuanet.com/content/2012-02/06/content_24649397.htm 


48 


http://dqnews.zjol.com.cn/dqnews/system/2010/08/17/012525402.shtml 


49 


http://cheshang.16888.eom/newsinfo/2011/1 1 15/1 41 264.html 
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Table 2. Cont. 



Website no. 


Website name of matchmaking fairs 


50 


http://bbs.heze.cc/thread-842865-1 -1 .html 


51 


http://www.8hy.Org/hyjy/hy6240/1 


52 


http://www.cdrb.com.cn/html/201 2-04/03/content_1 546492.htm 


53 


http://heilongjiang.dbw.cn/system/201 3/02/1 6/0545841 50.shtml 


54 


http://www.e0734.com/2012/0502/90707.html 


55 


http://sz.tznews.cn/tzwb/html/201 2-07/09/content_71 285.htm 


56 


http://www.xtrb.cn/epaper/ncwb/html/201 1-08/09/content_275667.htm 


57 


http://www.hukou365.com/cwbbs/forum/showtopic_tree.jsp7rootid = 194730 


58 


http://news.163.eom/1 0/0329/03/62TPSMMO0001 46BB.html 


59 


http://www.chinajilin.com.cn/content/2009-02/15/content_1495554.htm 


60 


http://sy.house.sina.com.cn/news/201 1-1 2-27/1 1 4483993.shtml 


61 


http://news.dayoo.com/guangzhou/201205/02/73437_23554932.htm 


62 


http://cq.qq.eom/a/20090824/000190.htm 


63 


http://news.qq.eom/a/201 1 1 228/000342.htm 


64 


http://www.efu.com.cn/data/201 1/201 1 -08-09/389729.shtml 


65 


http://ent.1 63.eom/1 2/1 203/1 3/8HQ885MN00032DGD.html 


66 


http://www.subaonet.com/html/society/2010426/3C95FFIB98JI5FC.html 


67 


http://wb.sznews.com/html/201 1-11 /07/content_1 81 2730.htm 


68 


http://heilongjiang.dbw.cn/system/2012/04/23/053819817.shtml 


69 


http://news.cnnb.com.cn/system/201 1/1 0/31 /007128083.shtml 


70 


http://www.people.com.en/GB/paper447/1 71 68/1 505082.html 


71 


http://news.timedg.com/201 2-04/1 6/content_9577975.htm 


72 


http://www.gddgart.com/artcenter/html3asp/town3ship/dq2012714_2357.asp 


73 


http://epaper.oeeee.eom/l/html/201 2-1 1/12/content_1 751 538.htm 


74 


http://news.hsw.cn/system/2010/06/28/050547974.shtml 


75 


http://fj.sina.com.en/news/s/201 2-08-24/071 86785.html 


76 


http://net.chinabyte.eom/1 64/1 221 01 64.shtml 


77 


http://epaper.hljnews.cn/shb/html/2008-05/26/content_199685.htm 


78 


http://www.estour.gov.cn/news/lvyouxinwen/2011/815/1181583840H7H40DE3ADE501J93BG9.shtml 


79 


http://www.0481 00.com.cn/news/bdxw/2009-04-20/61 7.html 


80 


http://www.estour.gov.cn/news/lvyouxinwen/2011/815/1181583840H7H40DE3ADE501J93BG9.shtml 


81 


http://www.estour.gov.cn/news/lvyouxinwen/2011/815/1181583840H7H40DE3ADE501J93BG9.shtml 


82 


http://zjnews.zjol.com.cn/05zjnews/system/2009/03/30/015385573.shtml 



doi:1 0.1 371 /journal.pone.0081 424.t002 



numbers of males and females would not differ over some extent. 
We narrow the range of r[ to [1,2] and solve this optimization 
problem 

min Q(n,)i), subject to 1 <)7<2,1 <n< 10 4 , ^elR, neZ. (14) 

Finally we obtain the estimation n= 15. 

Figure 3 shows the relationship between the experimental data 
and the analytical predictions of our model. The red curve and 
olive curve are obtained from (8). The parameters of red curve are 
v\= 1, n= 15; the parameters of olive curve are max(i7) = 365.67, 
n= 15. According to (7), when k\ is equal to the minimum 1, C 
takes the maximum value 225. The error bars of ordinate P of 
round dots represent the ranges of empirical data P in Table 1 . 
Because k\ is unknown and £ is undetermined, the bound for k\ in 
the z'th experiment is min(A;i ,) = E, <k\ ,<max(A:i iI ) = 



(kij+k2,i)/2, and the bound for corresponding is 
15 2 /max(^i ,)< q, < \5 2 /mm(k\ ,). Therefore, the ranges of 
abscissa C,- of round dots are relatively wide and the middle points 
lie in c, = (15 2 /max(/fc u ) + 15 2 /min(/c u ))/2. 

Figure 3 also shows when £ is relatively small and corresponding 
k\ is big, all empirical data are enclosed between two curves; when 
£ is relatively big and corresponding k\ is small, some empirical 
data are enclosed between the two curves, but other empirical data 
lie above the red curve and the trend of the empirical data is 
opposite to the analytical predictions. The possible reasons are: on 
the one hand, organizers of some matchmaking fairs select only a 
few participants meeting their requirements from a large number 
of applicants, so a participant is easier to find the right man or 
woman; on the other hand, when the number of participants is 
small in a matchmaking fair, they understand the difficulty of 
finding an ideal object so compromise to a goodish choice. The 
two reasons above cause that the fewer the participants are, the 
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Figure 3. The relationship between the experimental data and analytical predictions in the log-log plots. The red curve and the olive 
curve are obtained from (8), and the parameters of red curve are r\= 1, n= 15; The parameters of olive curve are ^ = 365.67, n= 15. The round dots 
represent the empirical data in Table 1. The £„ MV represents the maximum value 225 of 6 
doi:10.1371/journal.pone.0081424.g003 



higher the matching probability P is. Based on these effects, the 
deviation of experimental data from the model is acceptable. 

Conclusion 

We propose a model of the two-way selection system and 
provide its analytical solution. Under several conditions, the 
compact approximations are derived analytically and verified by 
the simulation results. In the model, the parameter n that denotes 
the number of characters directly determines the probability of the 
successful match - due to its importance, we propose a rough 
method to estimate its value by fitting the empirical data collected 



via the Internet and the result is M=15. Under some artificial 
assumptions, most of the experimental data fall into the range 
predicted by our model, so this model is helpful for understanding 
the dynamics mechanism of the real-world two-way selection 
systems, and provides a starting point for researching the nature of 
real-world two-way selection systems. We believe our model could 
enlighten readers in this rapidly developing field. 
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